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' Abstract 

^ ■ In this paper, we discuss the system of Friedman- Robertson- Walker (FRW) met- 

ric coupling with massive nonlinear dark spinors in detail, where the thermodynamic 

00 ' 

movement of spinors is also taken into account. The results show that, the nonlinear 
^2 ' potential of the spinor field can provide a tiny negative pressure, which resists the 

t^- ' Universe to become singular. The solution is oscillating in time and closed in space, 

which approximately takes the following form 

Qfiu = R 2 {1 — 5 cost) 2 diag(l, —1, — sin 2 r, — sin 2 r sin 2 0), 

with R = (1 ~ 2) x 10 12 light year, and 5 = 0.96 ~ 0.99. The present time is about 
t ~ 18°. 

PACS numbers: 98.80.-k, 98.80.Jk, 98.65.Dx, 95.35.+d 

Key Words: cosmological model, singularity- free, nonlinear 
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1 Introduction 

Combined with the fluid model, the general relativity usually gives the cosmological mod- 
els with singularity. After a series of singularity theorems, which is based on the conditions 
such as positive energy, closed trapped surface, have been proved by Penrose, Hawking 
and so on[l, 2, 3, 4], it seems that the singularity has been quite naturally accepted as an 
objective existence. However the intrinsic singularities are after all abnormal and incon- 
ceivable. The fluid model is just a simplified description of matter, so one is unnecessary 
to take it too serious if the intrinsic essences of the spacetime and matter are involved. 

*cmail: yqgu@fudan.edu.cn 
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The global solutions of Einstein's equation are always interesting[5]- [1 1] . The Lemaitre's 
Phoenix picture, a regular and oscillating cosmological model is reconsidered in [12, 13], 
it can give a more natural explanation for some famous cosmological puzzles such as the 
flatness and horizon problems. However this model is obtained by hypothesis rather than 
by solution of Einstein equation. Recently the coupling system of spinor field with space- 
time has been studied by several authors. The principal motivation of the works was 
to find out the regular solutions of the corresponding field equations[14]-[17] or explain 
the accelerating expansion of the present Universe [18]- [21]. These works indeed obtained 
singularity-free solutions, and showed that the nonlinear spinor field results in a rapid 
expansion of the Universe. 

In what follows, we also investigate the cosmological model with nonlinear spinor fields, 
which is derived from a general framework of field theory [22]- [25]. Different from the 
treatment of [15]- [21], in which the authors introduced only one spinor flied independent 
of spatial coordinates, and suffered from the energy momentum tensor incompatible with 
the high symmetry of the Friedman-Robertson- Walker(FRW) metric. In our point of 
view, one spinor field just describes one particle and the high symmetry is a result in 
mean sense. Then all problems can be easily solved, and the result is quite natural and 
reasonable. 



2 The Basic Equations 

A large number of empirical data shows that our Universe is highly homogeneous and 
isotropic. So the Universe can be described by FRW metric in mean sense, the corre- 
sponding line element is given by 

ds 2 = dr 2 - a 2 (r) ( dr * - + r 2 (d6 2 + sin 2 6d(p 2 )) , (2.1) 
V 1 — K r z J 

where K — 1, and —1 correspond to the closed, flat and open Universe respectively. 
However this form of FRW metric is not convenient for analysis, and the solution is 
difficult to expressed by elementary functions, so in this paper we adopt the conformal 
coordinate system[13]. The corresponding metric becomes 

9lUf = a 2 (t)diag (l, -1, -/ 2 (r), -/ 2 (r) sin 2 d) , (2.2) 

where dr = a(t)dt, 

{sin r if K — 1, 

r if K = 0, (2.3) 
sinh r if K — —I. 
The Lagrangian of the gravitational field is generally given by[26] 
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where R is the scalar curvature, A the cosmological factor. For the metric (2.2), we have 
the scalar curvature and Einstein tensor G^ u as follows 

a" + Ka . . 

R = 6—^ — , 2.5 

Goo = + (2.6) 

n" a' 2 

Gn = 2- - + K, (2.7) 

a ar 

G 22 = G n / 2 , G 3 3 = G 11 / 2 sin^, (2.8) 



where a' = |a, /' = J;/. 

Now we construct the Lagrangian of matter. Denote the Pauli matrices by 

'-*>-«: :)■(: ?)•(: '.))■ 

Define 4x4 Hermitian matrices by 

Instead of Dirac matrices 7 M , we adopt the above Hermitian matrices (2.10) for the conve- 
nience of calculation. Obviously, the most part of matter should be described by spinors. 
For some reasons, the following Lagrangian is the most natural candidate for matter fields, 
especially for dark matter. In flat spacetime the Lagrangian is given by 

C = E <i>t(°n% ~ ^l)ct> k + F{al p k , %), (2.11) 

k 

where we assign to the k-th dark spinor Sk, Hk > is constant mass, which takes one 
value for the same kind particles, 

«fe = 0fe" M 0fc ; h = 0t(3<i>ki % = <Ptl<Pk, (2.12) 

F is the nonlinear coupling term, which is small in value (\F\ <C Hk), but decisive for the 
structure of matter [24], [2 7] -[30] and the Universe. 

Considering that all the particle-like solutions of the nonlinear spinor equation have 
a mean diameter less than 10 2 Compton wave lengthes, and all \<frk\ decay exponentially 
with respect to the distance from the center [24, 29, 30], so we can neglect the nonlinear 
interaction terms among spinors, but keep the self coupling terms only. By the Pauli-Fierz 
identities = $ 2 + 7 2 [31] and the stability of the fields, F should take the following 

form 

F = Y. V ^ V k = V(%), (2.13) 

k 
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where the self-coupling potential V is a different iable, even and concave function satisfying 

V(0) = 0, V\x)x - V(x) > 0. (2.14) 

The covariant form of spinor equation in curved spacetime have studied by many 
authors [32]- [40]. The generalized form of (2.11) with diagonal metric is given by 

£m = E (pkfiWk) - Vk% + V k ) , (2.15) 
k 

in where 

fa a 1 a 2 a 3 \ . 
V a a aj aj sm # y 

The variation of (2.15) with respect to (j>l gives dynamic equation for Sfc[39, 40] 

^V M fe = - K)70fc, (Vfc) (2.17) 
where V M = (<9 M + T M ) , T M is the spinor connection 

r '= T i"* 9 ' °) ■ (2 - 18 » 

Coupling (2.4) with (2.15), we get the total Lagrangian of the Universe with dark 
spinors 

C = ~^{R - 2A) + E (ttlfid^k) - + V k ) . (2.19) 

k 

Generally the variation (2.19) with respect to g^ v gives the Einstein equation 

+ Ag^ + 8ttGT^ = 0. (2.20) 

In the spacetime with orthogonal coordinates, the energy momentum tensor T^ u is rela- 
tively simple [38, 40] 



E Q* {^UQ^ U + Q u ^n<Pk) + (KTfe - V h )sT) • (2.21) 



In (2.21) we substituted the dynamic equation (2.17) into C m . 

For the spinor at particle state[25, 41], we have the following classical approximation 



0+f/> fc - u^l-vl5(x-X k ), (2.22) 
iVV fc -> mnvPfa (2.23) 

where is the central coordinate of is its 4- vector speed. Then the classical 

approximation of (2.21) provides the energy momentum tensor of the ideal gas with a 
tiny negative pressure term 



= EK« + WkST) V 1 - v l 5 "^ - (2-24) 
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where Wk = Ir3(VI% ~ Vk)d 3 x > is the proper energy of contributed by Vk, which 
acts as negative pressure. The details of the classical approximation and local Lorentz 
transformation see [41, 42]. The numerical results show that < Wk <C mk, but it keep a 
spinor to be an independent particle[24, 29, 30]. In what follows we see how this wizardly 
tiny term resists the Universe from singularity. 



3 Solution to the Equation 

The complete Einstein equation (2.20) is an overdetermined system under the assumption 
of FRW metric, because the FRW metric only holds in mean sense. For an overdetermined 
system we should appeal to the variation principle. The action corresponding to (2.19) is 
given by 

I = fdt f Ca A dVt, (3.1) 



where dVt = f 2 sin Odrdddip is the angular volume independent of t. Variation (3.1) with 
respect to a, we get Euler equation as follows 

(3a" + 3Ka- 2Aa 3 ) Q + ^ f ^ (^y^) _ ^ % + 4 y fc ) a s dn = 0> (3 . 2) 

where Q is any angular volume large enough to contain sufficient particles, such that the 
statistical principle is valid. Substituting (2.17) into (3.2) we get 

a" + Ka - l\a 3 = ^ £ / (p k % + 3Vfa - AV k ) a 3 dn. (3.3) 
Define the proper mass of Sk by 



I . {Hklk + 3V£% - 4^4) d 3 x k , (3.4) 

J R 3 

where x k is the local Cartesian coordinates of the central coordinate system of S k - Ob- 
viously nik is a constant independent of t. Substituting it into (3.3) and making local 
Lorentz transformation [41], we get 

a" + Ka- ha* = ^ £ (3-5) 
The drifting speed Vk depends on a(t), solving the geodesic of Sk we have[43] 

v k = - r ^=, or yfl - v\ = — a = , (3.6) 
'a 2 + bi 



where b k is a constant determined by initial speed. Substituting (3.6) into (3.5), we get 

a» + Ka- 2 -Aa* = ^- £ ™ fcfl . (3.7) 
3 3U x k en J a 2 + b\ 
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Multiply (3.7) by a' and integrate it, again by (3.6) we have 



a +Ka - -Aa = — — ^ / = + Co- (3.J 

6 6U x k cn Jl - vt 



J k 

Now we examine the meaning of Cq. By (3.8), (2.5)-(2.8) and (2.21) we have 

C — n' 2 _l_ Kn 2 1 \n i — 8nG m k a 



„4t^0 1 mfcQ 



The equivalent pressure for stationary is given by [30] 

P k = ^(T ° - 7£) = ~(<f>ta id <l> k - » k % - 2Vl% + 3V k ), (3.10) 

where x k stands for the central coordinate system of S k . The numerical results show that 
for ground state 

/ P k d 3 x k <0. (3.11) 

In statistical sense, we have 

M k = J n ^Uig iV <p k )a 3 dQ 





■4 
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/i*> -A=f ((! + H 2 )(3P fc + + 2V£t* - 3H) - |% 2 (/i fc - V^)7 fc ) rf 3 x fc 



= /j*> [(3 + i$P fc + ft k % + 2V k % - 314 + v 2 k {Vl% - V k )} d 3 x k , 

where we also made local Lorentz transformation for the integrals of S k . Then we get 
T ° = jfr / n E fc (fe°zV o fc ) + ^7 ~ ^) a 3 ^ 



= ^Ex fee Q(M fc + / R3 (K7 fe -T4)Vl-^ 3 %) (3.12) 
= ^ Ex teS) J* 3 [(3 + + ^ k % + 3V k % - AV k ] d 3 x k . 

Substituting (3.4) and (3.12) into (3.9) we get 



C pa V 4^ / Pkd 3 x k < 0. (3.13) 

So Co is a negative quantity related to the nonlinear potential V. Of course (3.13) can 
not be used as definition in (3.8), because Gq and Tq are "overdetermined functions" with 
respect to FRW metric. 

The kinetic energy of S k is calculated by 

K k = ^— -m k , (3.14) 



Although the spinors are fermions, but if the mean distance among them is much larger 
than their mean radiuses, then they are all at particle states and the effects of exclusion 
principle actually vanish. So the kinetic energy of such spinors should satisfy the classical 
distribution, rather than Fermi-Dirac statistics. Considering the case of relativity, we 
have the following Maxwell-like distribution 



/(K)dK = W-^-exp (~) dK. (3.15) 



ir(kT) 3 ^\ kT 
Substituting (3.14) and (3.15) into (3.8), we get 

1 , „ 87rGa 



a> 2 + Ka 2 - Ua* = ^£ / (K + m k ) /(K)dK + Q 
3 6\l v~7 n J0 



x k en- 



E (lkT + m k )+C . (3.16) 



3^ x k en ^ 2 
In [43] we derived the relation T(a) as follows 

1 , m fhb 2 

-kT = 3.17 

2 5a(a + Va 2 + b 2 ) V ; 

where fa is mean mass of all particles, b is a constant determined by the initial data 



b 



5kT ( 5kT Q s 



a \ m \ 4m J 
Substituting (3.17) into (3.16) we finally get 



<1, (leV ~ 10 4 A). (3.18) 



_ / 352 \ 
a « + Ka 2 - -Aa 4 = 2R + a) + C , (3.19) 

where 

R = ^ E ™* (3-20) 

is the radius scale of the Universe, which is independent of a. 

In [44], the influences of the constants in (3.19) on the solution were analyzed in detail. 
The main results are that: (Rl). b ^ mainly influences the behavior as a —>■ 0, it has 
little influence on the present Universe. (R2). A mainly influences the behavior as a > R, 
it seems to be superfluous term without any effective use, so A = is the most natural 
choice. (R3). Eq.(3.19) has singularity- free solution as long as Cq < —^Rb. 

According to these results, we set A = b = for the following computation, the 
detailed computation will be given elsewhere. Then (3.19) becomes 

a' 2 + Ka 2 = 2Ra + C . (3.21) 
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Solving (3.21), we get the singularity-free solutions as follows 

fl(l-<J cost), (C = -R 2 (l-5 2 ), 0<6<1), if / = sinr, 
R{5 + |t 2 ), (C = -2R 2 5, 5>0), if / = r, (3.22) 

#(<5cosht- 1), (C = -R 2 (6 2 - 1), <S>1), if / = sinhr. 

Calculating the energy- momentum tensor by 8n GT^ U = —G^ v , we have 

3(1 - 25 cost + 5 2 ) m 5 2 -l „ „ 

r oo = p ^ n t t^-> = ^Vvl x 7f2' lf * = smr ' 3 ' 23 

87rG(l — o cost) z 87rG(l — o cost) 2 

3t 2 -2 
Too = 8vrG(5+|t^' Tll= 8vrG(5+|t^' if / = r ' ^ 

3(25cosht-l-5 2 ) 1-5 2 if f-sinhr (3 25) 

T °° " 87rG(l-5cosht) 2 ' Tn " 87rG(l-5cosht) 2 ' lf f ~ ^ ^ 

From the above equations we learn that, only for the closed case the energy Too is always 
positive, so only this solution is meaningful in physics. Moreover we find Tn < for all 
cases. Tn < is inconceivable for perfect fluid model, but from (2.21) and (2.24), we 
learn that it is normal for nonlinear spinor field. 



4 Sketch of the Universe 

In this section we associate the closed case in (3.22) with some observational data and 

give a sketch of the Universe. The detailed analysis will be given elsewhere. Define the 
cosmological time by r 

t= I a(t)dt = R(t-Ssmt). (4.1) 
Jo 

Then the age of the Universe T reads 

T = R(t a - 5 sint a ), (4.2) 

where t a stands for the age in angular coordinate system. The Hubble's parameter is 
defined by 

_ da_ _ a/_ _ 5 sin t a 

~ adr ~ a 2 ~ R(l - 5 cost a ) 2 ' 1 ' j 

The mass density p and the pressure P in usual sense are defined by 

= T o 3 l-25cost a + 5 2 
P 8ttCR 2 (l-5cost a ) 4 ' V ' J 

- 1 1 - r) 2 

1 87rGi? 2 (1 -5 costa) 4 ' ' 
Denote the critical density p c = ^qH 2 ~ 8 x 10 _30 (g/cm 3 ), we have 

p a 2 1 - 2£cost Q + 6 2 

^ = jr l+ ^= ( 5sint Q ) 2 • (4 - 6) 
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The equation of state is defined by 

w = P - = ^ . (4.7) 

p 3(1 -2<5cost a + 5 2 ) 1 ; 

By (4.7) we find w(t a ) < is an increasing function for < t a < tv, and 

1 + 5 

3(1 - b) 

So w of spinors is a function crossing the value —1. 

There are serval authoritative empirical data[45]-[53]. We take the data in [51] as basic 
parameters for computation. 

T = 14.2(Gyr), fi tot = 1.02 ~ 1.08, u = -1 0.2. (4.8) 

By computation we get 

t a = 12° ~ 22°, 5 = 0.96 ~ 0.99, R = (1.0 ~ 2.0) x 10 12 yr- (4-9) 

By t a we learn that the Universe is very young. The radius of the present Universe is about 
nR(l — 5 cos(t Q )) ~ 300Gyr. The period for a photon to travel a cycle in the Universe is 
just the life period of the Universe herself, so we can never observe the ghost of a galaxy 
at any time. However, by the functional digraphes of H(t), Q tot and w(t), we find they 
are rapidly varying functions of t of the past, then consequently , the observational data 
strongly depend on the distances of the galaxies. So the parameters provided here could 
only be used as rough scales, the more values should be obtained dynamically under the 
aid and constraint of specific cosmological models. 

From (3.22) we learn that, instead of the initial singularity there is a cradle of rebirth 
t = 0. At this time, the drifting speed of the particles or equivalently the cosmic tem- 
perature takes the maximum, and the volume of the Universe takes the minimum. This 
situation will result in violent collision among galaxies and particles. Such collisions will 
recover the vitality of the dead stars and combined particles. Just like the rebirth of the 
Phoenix, the Universe gets her rebirth via the destructive high temperature. 

Of course the above analysis are based on the idealized model of particles. What is 
important for this paper is that, it reveals the solutions of the Einstein's equation are 
sensitively related to the properties of the source. The interactional mechanism between 
the Universe and fundamental particles has been widely investigated by the group of the 
Center for Cosmoparticle Physics[54]-[66], and their methods and results will be important 
for further study on detailed interaction. 
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